OUTP 94-17 P 
HD-THEP-94-29 



Analytical Solutions of Exact 
Renormalization Group Equations 

N. Tetradis'^'Q and D.F. Litim'^'0 



^^Theoretical Physics 
University of Oxford 
1 Keble Road 
Oxford 0X1 3NP, U.K. 

'^Institut fiir Theoretische Physik 
Universitat Heidelberg 
Philosophenweg 16 
69120 Heidelberg, Germany 



Abstract 

We study exact renormalization group equations in the framework of the effec- 
tive average action. We present analytical solutions for the scale dependence of the 
potential in a variety of models. These solutions display a rich spectrum of physical 
behaviour such as fixed points governing the universal behaviour near second or- 
der phase transitions, critical exponents, first order transitions (some of which are 
radiatively induced) and tricritical behaviour. 
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1. Introduction 



The solution of an exact renormalization group equation is a particularly difficult 

task. The reason is that such an equation describes the scale dependence of an effective 
action, which is characterized by infinitely many couplings multiplying the invariants 
consistent with the symmetries of the model under consideration. As a result an exact 
renormalization group equation corresponds to infinitely many evolution equations for the 
couplings of the theory. The crucial step is developing efficient approximation schemes 
which can reduce the complexity of the problem while capturing the essential aspects of 
the physical system. Perturbative expansions have been used for proofs of perturbative 
renormalizability p, |ll|, while the powerful e-expansion ^ |12] has been employed for 



the study of fixed points governing second order phase transitions in three dimensions. 
In many cases evolution equations for truncated forms of the effective action have been 
solved through a combination of analytical and numerical methods (see the review 0] and 
references therein, |^ 0, [|T^-|^. Also, numerical solutions for the fixed point potential 
of three-dimensional scalar theories have been computed in ref. 



In this paper we present analytical solutions of equations which describe the evo- 
lution of the running effective potential. We are interested in describing the evolution 
not only near the possible fixed points, but also far away from them. This allows us to 
study the influence of the initial conditions (the parameters of the classical theory) on 
the renormalized theory. We can also investigate theories with first order phase transi- 
tions, which do not exhibit universal behaviour. The biggest part of this work concerns 
the 0(A^)-symmetric scalar theory in the large N limit, the spherical model ||23|. We 
obtain analytical solutions of the evolution equation which describes the dependence of 
the running effective potential Uk on an effective infrared cutoff k. This equation, and our 
solutions, are exact in the large N limit. In three dimensions, we show that for a certain 
range of classical parameters the evolution of the critical theory leads to an infrared at- 
tractive fixed point (the Wilson-Fisher fixed point). This results in a second order phase 
transition with universal behaviour. For a different parameter range the theory has a 
ffist order phase transition. A tricritical line separates the two regimes. In four dimen- 
sions we describe how the renormalized quartic coupling evolves to zero as the infrared 
cutoff is removed. Finally in two dimensions we explicitly demonstrate the absence of 
symmetry breaking and a phase transition for this model. The incorporation of physical 
phenomena of such a wide range is the most prominent and promising characteristic of 
our analytical solutions. We subsequently study an approximate evolution equation for 
the running effective potential in the context of the Abelian Higgs model in four dimen- 
sions. We demonstrate how its solution predicts a ffist order phase transition driven by 
radiative corrections (the Coleman- Weinberg mechanism). Parts of our solutions for the 
0(A^)-symmetric scalar theory have been obtained previously 0, |^. However, the start- 
ing point for our derivation is a recently proposed renormalization group equation |^, and 
our results demonstrate the complete agreement of its physical predictions with those of 
other renormalization group equations @, |^, |]. We also give a detailed discussion of the 
full range of the running effective potential (including its non-convex part), and consider 
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a variety of regions for the classical parameters of the theory. 

We work in the framework of the effective average action H, |l^, which has been used 
for the study of second and first order phase transitions for a variety of four- dimensional 
field theories at non-zero temperature |[T^-[|l^- Also, in ref. Jl^ an approximate solution 



iTr{(ri^) + i?,)-i|i?4. (1.1) 



of the evolution equation for the potential for the three-dimensional 0(A^)-symmetric 
scalar theory in the large limit has been given. The effective average action results 
from the integration of quantum fluctuations with characteristic momenta > k'^. It 
interpolates between the classical action S for k equal to the ultraviolet cutoff A of the 
theory (no integration of modes) and the effective action F for k = (all the modes are 
integrated). Its dependence on k is given by an exact renormalization group equation 
with the typical form (t = ln(fc/A)) 

Here F^ is the second functional derivative with respect to the fields, Rk is the effective 
infrared cutoff which prevents the integration of modes with <k'^, and the trace implies 
integration over all Fourier modes of the fields. We work with an approximation which 
neglects the effects of wavefunction renormalization. Therefore, only a classical kinetic 
term in the effective average action is kept, which takes, for an 0(A^)-symmetric scalar 
theory, the form 

Ffc = 1 rf"x{f/fe(0) + U^'ct>ad^r}, (1-2) 

and all invariants which involve more derivatives of the fields are neglected. With this 
approximation, eq. ( p..l|) can be turned into an evolution equation for the potential Uk- 
First we shall discuss the 0(A^)-symmetric scalar theory in the large A^ limit. In three 
dimensions the fixed point solution which governs the second order phase transition will 
be identified. We shall also show that, for a certain parameter range, the theory has a first 
order phase transition. The two regions in parameter space are separated by a tricritical 
line. We shall also study the same theory in four and two dimensions. Finally we shall 
discuss the Abelian Higgs model in four dimensions, for which the radiatively induced 
first order transition will be reproduced. 



2. The 0(A^)-symmetric scalar theory 

We first consider the evolution equation describing the dependence of the effective average 
potential Uk on the scale k in arbitrary dimensions ci, for an 0(A^)-symmetric scalar 
theory. As we have mentioned in the introduction, we neglect the effects of wavefunction 
renormalization. The evolution equation reads (t = ln(A;/A)) |^ 

dt 'Jo dt\P + U',{p)^ P + U',{p) + 2Ul{p)p 
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Here p 



1...N, and primes denote derivatives with respect to p. The variable 



X denotes momentum squared x = q , and 



(2.2) 



The inverse average propagator 



P{x) 



X 



incorporates an effective infrared cutoff k, such that modes with x 
propagate. The "cutoff" function /|(a;) is given by 



(2.3) 



-C k'^ do not 



/fc(x) = exp<^ -2a 



X 



(2.4) 



and can be made sharper or smoother by an appropriate choice of the two free param- 



eters a, b. Up to effects from the wavefunction renormahzation eq. ( |2 . 1| ) is an exact 
non-perturbative evolution equation It is easy to recognize the first term in the r.h.s. 
of eq. ( |2 . 1| ) as the contribution of the — 1 Goldstone modes (f/(, vanishes at the mini- 
mum). The second term is related to the radial mode. After performing the momentum 
integration the evolution equation ( p . 1|) becomes a partial differential equation for Uk with 
independent variables p and t. The effective average potential interpolates between the 
classical potential V for k = A (with A the ultraviolet cutoff) and the effective potential U 
for /c = [§]. As a result the solution of eq. (|2.1|) with the initial condition Ua{p) = V{p) 
uniquely determines, for /c — 0, the effective IPI vertices at zero momentum for the 
renormalized theory. In order to write eq. ( |2.1|) in a scale invariant form it is convenient 
to define the variables 



p =e~'^p 

Uk{p) =k-''Uk{p). 



(2.5) 



In terms of these eq. ( |2.1|) can be written as (for the details see refs. p], |T5|, |T6|, IT 

du' 



dt 



- 2u' +{d- 2)pu" + Vd{N - l)u"L'i{u') + Vdi'iu" + 2pu"')Ll{u' + 2pu"). (2.6) 



Primes on u denote derivatives with respect to p and we omit the subscript k in from 
now on. The functions Lf[w) are given by 



Lf(w) = -2(2a)— / dy y—e'^ 
Jo 




(2.7) 



(where we have used the variable y = 2a {x/k"^)^). They introduce "threshold" behaviour 
in the evolution equation, which results in the decoupling of heavy modes. They approach 
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a constant for vanishing argument, and vanish for large arguments. For general a and b 
they do not have a simple analytical form. However, in the sharp cutoff limit b ^ oo they 
are given by the simple expression 

LfW= (2.8) 

The resulting evolution equation (with an appropriate rescaling of the fields) is identical 



to the one introduced in ref. B and studied in refs. B, 0, |2^. Also, in the large N 



limit (see below) it is the same as the evolution equation which results from the exact 



renormalization group equation of ref. [Q, for which solutions are given in ref. p^ |. 
Our derivation establishes the agreement of the physical predictions resulting from the 
formalism of the effective average action with those of previous renormalization group 
formulations. We should point out that the use of a sharp cutoff renders the effective 
average action Tk non-local. This complication, however, does not affect our analysis, 
as we have neglected the effect of the higher derivative terms on the evolution equation 
for the potential. For studies which take into account the effects of the wavefunction 
renormalization Jl^ the use of a smooth cutoff is necessary. 



Even with the simple form of eq. (|2.8| ) for the "threshold" functions, the evolution 
equation ( p.6|) remains a non-linear partial differential equation which seems very difficult 
to solve exactly. The problems arise from the contribution from the radial mode to the 
r.h.s. of eq. ( |2.6|) (the last term). An enormous simplification is achieved, however, if this 
contribution is neglected. This is justified in the large limit, in which the contributions 
of the Goldstone modes dominate the evolution, and the contribution of the radial mode 
becomes a sub leading effect. Another important simplification results from the fact that 
the anomalous dimension of the field is zero to leading order in 1/A^, in four and three 
dimensions which are of most interest As a result, our approximation of neglecting 
the wavefunction renormalization is justified in this limit and our solution of eq. ( |2.6|) 
becomes exact. We should point out that the contributions of the radial mode cannot 
be neglected even in the large limit in some cases. For certain forms of the classical 
potential, the approach to convexity for the effective potential depends crucially on the 
radial mode. An example of this behaviour will be given in the following. For A^ 1 the 
simplified evolution equation reads 

du' , , ^dv! NC du' ^ , ^ 

with C = 2vd- It is first order in both independent variables and can be solved with the 
method of characteristics. 

a) Three dimensions 

We are interested in the behaviour of the theory in three dimensions, where a non-trivial 
fixed point structure arises. The most general solution of the partial differential equation 



■^See ref. Q for ^ numerical solution of eq. (2.6) 
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(|2.9D for (i = 3 is given by the relations 



P 



NC 



+ NC arctan 



F u'e 



l2t 



NC 



-NC In 
2 



1 + 



F(u't 



2t 



for m' > 
for u' < 0, 



with 



C = 2v3 



1 



47r2 



(2.10) 
(2.11) 

(2.12) 



The function F is undetermined until initial conditions are specified. For t = {k = A), 
Uk coincides with the classical potential V. The initial condition, therefore, reads 



u'{p,t = 0) = A-'V'{p). 
This uniquely specifies F and we obtain 



(2.13) 



p NC 



'u 
G (u'e"' 



+ NC arctan 
NC 



+ NC arctan 



for m' > 



(2.14) 



P 



NC 



2 U + v^, 



NC 



\/— w'e* \J—u'e* 2 



iVCln 



1 - V^^/e* 
1 + y^e* 



for m' < 0, 



(2.15) 



with the function G determined by inverting eq. (|2.13|) and solving for p in terms of u' 

G{u') = p{u%=o. (2.16) 

We are also interested in the scale invariant (independent of t) solutions of eq. (|2.9|) . 
These solutions correspond to the possible fixed points of the evolution and are obtained 
by setting the first term in the l.h.s. of eq. ( |2.9|) equal to zero. We can easily identify the 
"trivial" solution = 0, which is independent of the number of dimensions. For d = ?> 
the only other solution which remains finite for finite p is given by the relations 



P 



NC 



P 



+ NC arctan 



7r 



NC 



NC 



1 



-NC In 



1- 







for u^>0 
for u[ < 0. 



(2.17) 
(2.18) 



This solution corresponds to the Wilson-Fisher fixed point. A large part of the above 
results has been obtained in the past |]14|, 0. However, we have explicitly presented here 
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the solution for the non-convex part of the potential. We are also mainly interested 
in studying the various phase transitions which result for different ranges of classical 
parameters of the theory. 

I) Classical (p^ theory: Let us first consider a theory with a quartic classical potential. 
The initial condition can be written as 

«'(p,t = 0) = Aa(p-ka), (2.19) 

with 

= ^. A» = ^ (2.20) 

the rescaled (dimensionless) minimum of the potential and quartic coupling respectively. 
The function G in eqs. ( |2.14|) , ( |2.15|) is now given by 

Gix) = K^ + ^. (2.21) 
The typical form of the effective average potential Uk{p) at different scales k, as given by 



eqs. (|2.14| ), ( |2.15|) , is presented in fig. 1. The theory at the ultraviolet cutoff is defined 
in the regime with spontaneous symmetry breaking, with the minimum of the potential 
at poA = ^aA 7^ 0. At lower scales k the minimum of the potential moves continuously 
closer to zero, with no secondary minimum ever developing. We expect a second order 
phase transition (in dependence to k,\) for the renormalized theory at k = 0. There is a 
critical value for the minimum of the classical potential 

= i^cr = NC, (2.22) 



for which the scale invariant (fixed point) solution of eqs. (|2.17|) , (|2.18| ) is approached 



in the limit t — oo {k 0). Eqs. ( |2.17| ), ( |2.18|) describe a potential u which has a 
minimum at a constant value 

K{k) =K^ = NC. (2.23) 

This leads to a potential Uk{p) with a minimum at po(^) = kn^, —>■ for k —>■ 0, which 
corresponds to the phase transition between the spontaneously broken and the symmetric 
phase. (The values for Hcr and ft* coincide, but this is accidental.) For the second and 
third p-derivative of u at the minimum A = u"{k), a = u"'{k) we find 

m =A. = (2.24) 

and similar fixed point values for the higher derivarives of u. For 1 <^ p/NC ^ (7r/2)e^* 
the rescaled potential u has the form 
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Notice that the region of vahdity of eq. ( |2.26| ) extends to infinite p for t — oo. From 

eq. (p.26|) with t — > — oo(A; 0) we obtain for the effective potential at the phase 
transition 

u 2 ^^3 



Through eqs. ( |2.14| ), (|2.15| ) we can also study solutions which deviate slightly from 



the scale invariant one. For this purpose we define a classical potential with a minimum 

KA = Kcr + ^^A, (2.28) 

with \5k\\ <^ NC. A typical solution is depicted in fig. 2. The potential starts its 
evolution with its classical form at /c = A. It subsequently evolves towards the scale 
invariant solution given by eqs. ( |2.17|) , ( |2.18|) . It stays very close this solution for a long 



"time" t, which can be rendered arbitrarily long for sufficiently small |(5ka|. Eventually it 
deviates towards the phase with spontaneous symmetry breaking or the symmetric one. 
In fig. 2 (5ka was chosen negative, so that the evolution for A; — > leads to the symmetric 
phase. The curvature of the potential at the origin becomes positive and u'(0) diverges so 
that the renormalized mass term U'f.{Q) = k'^u'{0) reaches a constant value. The evolution 
of the minimum of the potential is given by 

K{k) = + ^KAC"*. (2.29) 

Also for A we find 

m = (2.30) 

1 + - 1) e* 



Eq. ( |2.29| ) indicates that the minimum of u stays close to the fixed point value given by 
eq. ( |2.23| ), for a very long "time" \t\ < — In |5ka|- For \t\ > — In \ 5kx\ it deviates from the 
fixed point, either towards the phase with spontaneous symmetry breaking (for 6k\ > 0), 
or the symmetric one (for 5ka < 0). Eq. ( 2.30 ) implies an attractive fixed point for A, 



with a value given by eq. ( |2.24| ). Similarly the higher derivatives are attracted to their 
fixed point values. The full phase diagram corresponds to a second order phase transition. 
For 6k,a > the system ends up in the phase with spontaneous symmetry breaking, with 

Po = lini Po{k) = lim kn{k) = 6k,aA. (2-31) 

In this phase the renormalized quartic coupling approaches zero linearly with k 

Xr = lim kX{k) = lim kX^ = 0. (2.32) 

The fluctuations of the Goldstone bosons lead to an infrared free theory in the phase with 
spontaneous symmetry breaking. For 6k,\ < 0, K,{k) h GCOIXIGS zero 8bt cL SCclIg 
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and the system ends up in the symmetric regime (po = 0). From eq. (|2.14D, in the hmit 



t — > — oo, with u',u",u"' oo, so that w'e^* ~ jf^fi^Ap, u"e^ 



U{p) = Uoip) 



TC 



2 y 

Nc) 



1, we find 



(2.34) 



Notice how every reference to the classical theory has disappeared in the above expression. 
The effective potential of the critical theory is determined uniquely in terms of Ska, which 
measures the distance from the phase transition. The above results are in exact agreement 
with refs. |14|, |l^, In particular the values for the critical exponents /?, z/, 

describing the behaviour of the system very close to the phase transition, correspond to 
the large N limit of the model (for details see section 10 of ref. [p!7| ) 

d (\n^ 



(2.35) 



P= lim ; ( = 0.5 

<5ka^o- (i(ln I^kaI) (i(ln|(5KA|) 



We also point out that, for a theory with spontaneous symmetry breaking , we can use 
eq. ( |2.15|) in order to study the "inner" part of the potential. In particular, for p = and 
t — > — oo eq. (|2.15|) predicts a potential u which asymptotically behaves as 



lim u'(0) 



-1. 



This leads to an effective average potential Uk which becomes convex with 

limf/^(0) = -A;^ 

fe— >o 



(2.36) 



(2.37) 



in agreement with the detailed study of ref. 

II) Classical (j)^ theory: As a second example we consider a theory defined through a 
classical potential with a p^ (0^) term 



u\p, t = 0) = Aa(p - /tA) + ^(p - I^Kf, 



(2.381 



where ka, Aa are defined in eq. (|2.20| ) and the coupling a a is dimensionless in c? = 3. The 
function G in eqs. ( |2.14D , ( |2.15| ) is now given by 



G{x) =ka + 



-Aa + \/Aa + So-ax 

f^A 



G{x) =ka + 



X\ + 2oax 



for u" > 



for u" <0. 



(2.39) 



(2.40) 



We distinguish two regions in parameter space which result in two different types of 
behaviour for the theory: 
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(a) For < '^^a/cta the classical potential has only one minimum at poA = kaA. Near 
this minimum the initial condition of eq. ( p.38|) is very well approximated by eq. ( p.l9|) . 
As a result, for ka near the critical value of eq. (|2.22D , the critical theory has exactly the 
same behaviour as for a quartic classical potential. The running potential first approaches 



the fixed point solution of eqs. (|2.171 ), (|2.18| ) (notice that < 2A^/crj,), and subsequently 



evolves towards the phase with spontaneous symmetry breaking or the symmetric one. 
The behaviour of the critical theory for k = is determined only by the distance from the 
phase transition (as measured by 5/€a), without any memory of the details of the classical 
theory. This is a manifestation of universality, typical of second order phase transitions. 
For ka > 2Aa/o"a the classical potential has two minima, one at the origin and one at 
PoA = kaA. The minimum at the origin is less deep for ka < SAa/cta. Again, for ka near 
Kcr the scale invariant solution is approached. We demonstrate this type of behaviour in 
fig. 3 for a theory with ka = Kcr- The classical potential includes a (p^ term and has a (less 
deep) second minimum at the origin, but the universal scale invariant solution is again 
approached for k ^ 0. Small deviations of ka from the critical value result in universal 
behaviour for the renormalized theory. 

(b) The minimum of the classical potential at the origin is deeper for k,\ > 3X\/(T\. 
An example of the evolution of the effective average potential Uk{p) for such a theory is 
given in fig. 4. The minimum of the potential at non-zero p moves towards the origin for 
decreasing scale k. In the same time the positive curvature at the origin decreases. The 
combined effect is (very crudely) similar to the whole potential being shifted to the left of 
the graph. As a result the minimum at the origin becomes shallower. For a certain range 
of the parameter space (for small enough k\, such as chosen for fig. 4) the minimum away 
from the origin becomes the absolute minimum of the potential at some point during the 
evolution. This results in a discontinuity in the running order parameter. Finally the 
absolute minimum of the potential settles down at some non-zero po. For even larger ka 
the minimum at the origin is deep enough for the evolution to stop while this minimum 
is still the absolute minimum of the potential. When the minimum of the renormalized 
potential po (which is obtained at the end of the evolution) is considered as a function of 
K\, a discontinuity is observed in the function Po(ka). This indicates a first order phase 
transition. Unfortunately, an exact quantitative determination of the region in parameter 
space which results in first order transitions is not possible within the approximations we 
have used. The reason for this is the omission of the "threshold" function for the radial 
mode which includes the term 2pM". As a result our approximation is not adequate for 
dealing with the shape of the barrier in the limit k —>■ 0, where the theshold function for 
the radial mode becomes important. Also the approach to convexity cannot be reliably 
discussed (in contrast to the case of a classical 0"^ potential). If the shape of the barrier 
cannot be reliably determined the relative depth of the two minima cannot be calculated, 
and our discussion is valid only at the qualitative level. 

However, more information can be extracted from our results. As long as we concen- 
trate on regions of the potential away from the top of the barrier the solution given by 
eqs. (|2.14|) , ( |2.39|) , (|2.4CI|) is reliable. This means that we can study the potential around 
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its two minima. We are interested in the limit t — oo (/c — > 0), with U' = w'e^*, p = pe* 
approaching finite values. The form of the potential near the minimum away from the 
origm IS determined by eqs. (|2J|), ( ^^391) . We find 



p -Aa + \ Xl + 2a [U' I \ 

E^-KA + NC = ^ ^ + iVCi/^arctan ^ . (2.41) 






(7a 

The minimum po (where U'{po) = 0) is located at po = (ka — A^C)A = 6haA. This requires 
6k,a > 0. Eqs. (p.l4| ), ( p.40|) describe the form of the potential around the minimum at 
the origin. Similarly as above we find 

^-nA + NC= ^2 ' ^"^^'^ + ^^\/^ -^^t-- I ^Tf I • (2-42) 



In the parameter range ka — 2Aa/o"a, 2Aa/o"a ^ NC the above solution reproduces the 
classical potential, with a large positive curvature ?7'(0)/A^ at the origin. This is due 
to the fact that the fiuctuations which renormalize the potential around the origin are 
massive, with their masses acting as an effective infrared cutoff. For the above parameter 
range these masses are of the order of the ultraviolet cutoff A and no renormalization 
of the potential takes place. This is in contrast with the form of the potential near the 
minimum po away from the origin. The presence of the Goldstone modes in this region 
always results in strong renormalization. There is a range of parameters for which the 
curvature at the origin becomes zero. It is given by the relation 

KA = NC + 2AA/fTA. (2.43) 

For this range the minimum at the origin disappears and the potential has only one 
minimum at po = SkaA = {ka — NC)A = 2Xa/(Ta- The above condition does not 
determine precisely the first order phase transition, as this takes place when the two 
minima are degenerate, and not when the minimum at the origin disappears. However, 
it provides a good estimate of its location. The discontinuity in the order parameter is 
expected to be 0{6ka)- Weakly first order transitions are obtained for Aa 0. We should 
emphasize that eq. ( |2.42|) is not valid for arbitrarily small f/'/A^. This would correspond 
to a range of the potential near the top of the disappearing barrier, where we know that 
our approximation fails. This is another reason why eq. ( |2.43| ) is only indicative of the 
location of the first order phase transition. 

We have identified two critical surfaces in parameter space. We saw in (a) that the 
surface ka = NC corresponds to second order phase transitions. Also in (b) we argued 
that the surface ka = NC + 2AA/crA corresponds to first order transitions. As a result we 
expect tricritical behaviour to characterize their intersection, which is given by the line 
Ka = NC, Aa = 0. This is confirmed if we approach this line close to the critical surface 
Ka = NC. More specifically we consider a theory with < —Ska = —ha+NC <C NC and 
Aa <^ 1/NC. For this choice of parameters the renormalized theory is in the symmetric 
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phase very close to the second order phase transition. The form of the potential is given 
by eq. (|2^ with U' / < 1 



A 




(2.44) 



For \5kx\/NC <^ \\NC the potential has the universal form of eq. ( |2.34| ). The initial 
point of the evolution is sufficiently close to the critical surface for the flows to approach the 
Wilson-Fisher critical point before deviating towards the symmetric phase. The critical 
exponent v takes the large N value u = 1 according to eq. (|2.35|) . In the opposite limit 
\6k\\/NC ^ XaNC the potential near the origin is given by 



U{p) = Aa 



(2.45) 



and the exponent u takes its mean field value u = 0.5. A continuous transition from one 
type of behaviour to the other (a crossover curve) connects the two parameter regions. 
Clearly, the line = NC, Aa = gives tricritical behaviour with mean field exponents, 
in agreement with the analysis of ref. |26|. R 



b) Four dimensions 

In the previous subsection we explored the non-trivial fixed point structure of the three- 
dimensional scalar theory. We now look for similar structure in four dimensions. The 
most general solution of the partial differential equation ( |2.9|) for d = 4 is given by the 
relation 

p NCI NC/l+u'' 



u' 2 u' 2 V \u 

where now 



=^K)' (2.46) 



C = 2v, = — ^. (2.47) 



We consider a quartic classical potential according to eq. (|2.19|) , with 

f^A = Aa = Aa. (2.48) 

This uniquely specifies F and we obtain 



p NCI NC^ fl + u'\ 1 / NC\ 1 NC ^ fl + u'e 



2 u' ' 2 ~"\ \u'\ J Aa ' V"'' 2 7 w'e^* 2 \ \u'\e^^ J' 

(2.49) 

There is again a critical value for the minimum of the classical potential 

NC 

K\ = Kcr = — r— (2.50) 



*The Bardeen-Moshe-Bander phenomenon [E7| was not considered in our discussion. 
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which separates two possible phases. The scale invariant solution is obtained for the above 
value of fCA in the limit t — > — oo. The resulting solution is the "trivial" one 

< = 0. (2.51) 

As a result, no interesting universal (independent of the classical parameters) structure 
can be obtained for this theory in four dimensions. 

We can also investigate solutions which deviate slightly from the critical one, by con- 
sidering a classical potential with a minimum given by eq. (|2.28|) . The renormalized 



potential in the limit t — > — oo, with m' — oo and finite U' = w'e^*, is given by the relation 



gH-^A.|;i„(l±£]=A.(|j-fc.) (2.52) 



for U' > 0. For theories with characteristic mass scale much smaller than the ultraviolet 
cutoff (such as in the symmetric phase in the vicinity of the phase transition) we have 
f/'/A^ — >■ 0. This gives for the renormalized quartic coupling 



The same happens near the minimum of the potential in the phase with spontaneous 
symmetry breaking, where no renormalized mass can be generated for the radial mode. 
All this is in agreement with the arguments for the "triviality" of the scalar theory in the 



limit that the ultraviolet cutoff is removed (see for example ref. and refs. therein). 



c) Two dimensions 

We finally turn to two dimensions, where the most general solution of eq. ( |2.9|) is given 
by 

with 

C = 2v2 = — (2.55) 

4:71 

and p dimensionless. We consider a quartic classical potential according to eq. ( |2.19| ), 
with 

i^A = PoA, Aa = (2.56) 

This leads to 

«'e- = A.(p-K.) + ^A.l„(|i±gj). (2.57) 

For any choice of classical parameters the minimum of the effective average potential runs 
to zero at a scale 

t. = (2.58) 



12 



At t = the renormalized potential has a minimum at zero and is given by the expression 

^ = Aa (p - K^) - In [jfw) ■ (2-59) 

The renormahzed theory cannot be defined in the phase with spontaneous symmetry 
breaking and there is no phase transition, in agreement with the Mermin- Wagner theorem 
21. 



3. The Abelian Higgs model in four dimensions: 

We now turn to gauge theories, for which exact renormahzation equations have also been 
obtained [^. As an example we discuss the Abelian Higgs model with real scalars, in 
four dimensions. The evolution equation can be written in the form (for the details see 
ref. m) 



^ = -2u' + 2pu" + {N- l)viu"L\{u') + v^{?,u" + 2pu"')L{{u' + 2pu") + Qv^e^L{{2e^p), 

(3.1) 



with V4 given by eq. (|2.2D. We have again neglected the small wavefunction renormal- 



ization effects for the scalar field. We recognize the contributions of the Goldstone and 
radial modes. The last term in eq. (|3.1| ) is the contribution of the gauge field. It involves 
the gauge coupling e^, whose evolution can be computed independently [§, |19]. Since 



the resulting running for is only logarithmic in d = 4, it is a good approximation to 
neglect it in the following and use a constant e^. The "threshold" functions Lf are given 



by eq. (|2.8|) . The contribution of the radial mode introduces higher derivatives in the 
evolution equation, making an explicit solution impossible. We shall again neglect this 
contribution, as in the first part of the paper. We should point out, however, that the 
resulting approximate evolution equation does not become exact any more in the large 
A^ limit. The purpose of this section is simply to demonstrate that the correct physical 
behaviour is incorporated in the full evolution equation, even though the approximate 
equation that we are using does not permit quantitative accuracy. After the omission 
of the contribution of the radial mode, the resulting partial differential equation is first 
order and can be solved with the method of characteristics. We have not managed to 
obtain an analytical solution in closed form, even though a numerical solution is possible. 
For this reason we make an additional approximation which is not crucial for the physical 
behaviour that we are interested in (see below). We set Lf{u') = Lf{0) = —2 in the con- 
tributions of the Goldstone modes, while maintaining the full "threshold" function in the 
contribution of the gauge field. As a result we cannot observe the decoupling of the scalar 
modes or the approach to convexity for the effective potential. However, we preserve the 
full effect of the gauge field on the form of the potential. We thus finally arrive to the 
following evolution equation 

du' r.~du' ^du' De^ „ , ^ /„ 
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with 

A^- 1 



B =2{N - l)v4 



D =l2vi = (3.3) 



The most general solution of eq. ( |3.2| ) is given by 
u' De^ 1 De^ , f2e^p + l 



2p-B 2{Be^ + l)2p- B 2{Be^ + 1)^ \\2p - B\ 

The function F is determined through the initial condition for the potential. Assuming a 
quartic classical potential given by eqs. ( |2.19|) , (|2.48|) we find 

^, , ^ x + B ^ 1 De^ 1 De^ , f e\x + B) + 1\ 

In fig. 5 we plot the potential which results from eqs. (|3.4| ), ( |3.5|) for a certain choice of the 
parameters of the theory. Initially the effective average potential has only one minimum 
at a non-zero value of p. As k is lowered a second minimum appears around zero, which 
eventually becomes the absolute minimum of the potential. The discontinuity in the 
expectation value signals the presence of a first order phase transition. The development 
of the minimum around zero is caused by the logarithmic terms in eqs. ( |3.4| ), ( |3.5| ). 
The situation is typical of a Coleman- Weinberg phase transition triggered by radiative 
corrections The effective potential U = Uq can be calculated from eqs. ( |3.4| ), ( |3.5| ) 
in the limit t — oo, with u', p oo, so that w'e^* ~ 1, pe^* ~ 1. We find 



A2 



p B De^ 



A2 \ 2 2AA(5e2 + l) 



(5e2 + l)2A2 \^e2(2^ + 5) + l 

(s'.e) 

Without the logarithmic term the phase transition in dependence to would have been 
second order. The presence of the last term in eq. ( |3.6| ) results in the development 
of a barrier near a secondary minimum at the origin. This leads to a weakly first order 
transition, with a discontinuity for the expectation value much smaller than the minimum 
of the classical potential. (For a detailed discussion of the Coleman- Weinberg transition 
using the full evolution equation see ref. |jl9|.) The effective potential of eq. ( |3.6|) is not 
convex. As we have mentioned already, the reason for this is the approximation of the 
"threshold" function by a constant in the evolution equation. Starting from eq. ( |3.6|) 
we can derive a relation between the different mass scales in this model. We denote the 
scalar field mass at the minimum po of the potential by m? = 2U"{po)po, the gauge field 
mass by = 2e'^po, and the scalar mass around the origin by mj. = U'{0). We find 

M2 " 87r2 (5e2 + l){Be^ + 1 + / K^) ' ^ ' ' 



The well-known result of Coleman and Weinberg is obtained for M2/A -C 1, if we 
neglect the scalar contributions by setting 5 = 0, and consider a potential with zero 
curvature at the origin by setting ml = 0. 
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4. Conclusions 



We presented analytical solutions of the exact renormalization group equation for the 
effective average action which are not restricted in the vicinity of a possible infrared 
attractive fixed point. This allowed us to investigate universal and non-universal aspects 
of phase transitions for a variety of models. We neglected the effects of wavefunction 
renormalization and approximated the action by the potential and a standard kinetic 
term. We solved the evolution equation for the potential as a function of the field and the 
running scale k. We presented analytical solutions for the 0(A^)-symmetric scalar theory 
in the large N limit in three, four and two dimensions. The omission of the effects of 
wavefunction renormalization is justified in four and three dimensions by the vanishing 
of the anomalous dimension to leading order in Parts of our results have been 

obtained in the past [|1^, 0]. However, the emphasis of our work lies in the detailed study 
of the influence of the classical parameters of the theory on the possible phase transitions. 
We also gave a crude treatment of the Abelian Higgs model in four dimensions. The 
enhanced complexity of the evolution equation for this last model made necessary the use 
of additional approximations, which, however, do not affect the qualitative conclusions. 

a) For the 0(A^)-symmetric scalar theory in the large N limit in three dimensions we 
distinguish two types of behaviour: I) For a classical 0^ potential given by eq. ( p.l9|) the 
renormalized theory has a second order phase transition in dependence on ka. The uni- 
versal behaviour near the transition is governed by the Wilson-Fisher fixed point and can 
be parametrized by critical exponents. II) For a classical (j)^ potential given by eq. ( |2.38|) 
there is a parameter range for which the renormalized theory has a second order phase 
transition in dependence on ka, with universal critical behaviour. For another parameter 
range the theory has a first order phase transition. The two regions are separated by a 
tricritical line (at Aa = 0) which displays tricritical behaviour with mean field exponents. 

b) For the same theory in four dimensions the solution reproduces the "triviality" of the 
critical theory and the vanishing of the renormalized quartic coupling. 

c) In two dimensions the solution demonstrates that the renormalized theory cannot be 
defined in the phase with spontaneous symmetry breaking and there is no phase transi- 
tion, in agreement with the Mermin- Wagner theorem. 

d) For the Abelian Higgs model in four dimensions we reproduce the Coleman- Weinberg 
first order phase transition which is triggered by radiative corrections. 

Our results on the universal behaviour of the three-dimensional scalar theory and on 
the four-dimensional Abelian Higgs model are in agreement with refs. [jl6|, |l^ p!9| . 



where the evolution equations for the effective average action were studied with other 
methods. Thus they provide an additional argument for the validity and applicability of 
the method of the effective average action in a wide range of problems. The most impor- 
tant aspect of our solutions, however, is that they are fully analytical and not restricted 
in the vicinity of a possible infrared fixed point. They provide a concise, transparent pic- 
ture of universal and non-universal behaviour, at all values of the field and the effective 
infrared cutoff, for a variety of physical systems. 
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Figure captions 



Fig. 1 The effective average potential Uk{p) at different scales for a scalar model in the 



large limit. The classical potential is given by eq. (|2.19|) with ka = 1.2NC and 
Aa = 0.3/A^C {d = 3). 

Fig. 2 The derivative u'{p) of the rescaled potential at different scales for a scalar model 
in the large N limit. The classical potential is given by eq. ( |2.19|) with slightly 
smaller than NC and Aa = 0.3/ NC {d = 3). 

Fig. 3 The derivative u'{p) of the rescaled potential at different scales for a scalar model 
in the large N limit. The classical potential is given by eq. ( |2.38|) with k\ = NC, 
Aa = 0.3/A^C and cta = OM/{NC)'^ {d = 3). 

Fig. 4 The effective average potential Uk{p) at different scales for a scalar model in the 
large N limit. The classical potential is given by eq. ( |2.38|) with ka = 1.2NC, 
Aa = 0.3/iVC and cta = l/{NCy {d = 3). 

Fig. 5 The approximate effective average potential Uk{p) at different scales for an Abelian 
Higgs model with = 2 and = 11 real scalars. The classical potential is given by 
eq. ( ^TTqI) with ka = 1 and Aa = 0.01 (rf = 4). 
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